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Combined D0 and CDF Upper Limits on Standard-Model Higgs-Boson Production 

The TEVNPH Working Group* 
for the D0 and CDF Collaborations 

We combine results of CDF and D0 searches for a Standard-Model Higgs boson (H) in data 
from pp collisions at the Fermilab Tevatron with y/s = 1.96 TeV. With 260-950 pb" 1 collected at 
VO _1 D0 ,and 360-1000 pb _1 collected at CDF, the 95% CL upper limits are a factor of 10.4(3.8) higher 

than the expected cross section for ran =115(160) GeV/c 2 . This result extends significantly the 
. individual limits of each experiment. 

(N ; 
O ■ 

<D ■ I- INTRODUCTION 

Q: 

qq ' Because the mechanism for electroweak symmetry breaking has yet to be confirmed, the search for a Standard- 

T — I , Model (SM) Higgs boson represents a significant issue in Fermilab's Tevatron physics program. Both D0 and CDF 
have recently reported searches for Higgs bosons that combined different final states and production modes |l|, |2j. 

In this note, we combine results of all such searches from CDF and D0 for pp collisions at sqrts — 1.96 TeV. These 
searches are for SM Higgs bosons produced in association with vector bosons (pp — > W/ZH —> £vbb/vvbbj '£ + £~bb or 
pp — > WH — * WW + W~) or singly through gluon-gluon fusion (pp — > H — > IF+TF"). The results are for data 
corresponding to integrated luminosities ranging from 360-1000 pb _1 at CDF and 260-950 pb _1 at D0 . The searches 
' are separated into sixteen final states, referred to as analyses in the following. To simplify combination of signals, the 

— , , analyses are separated into sixteen mutually exclusive final states. Selection procedures for each analysis are detailed 

\& • in each of the experiment's reports [l], H|, and briefly described below. 

O 

5_> 1 II. ACCEPTANCE, BACKGROUNDS, AND LUMINOSITY 

Oh" 

\ Selections are similar for the corresponding CDF and D0 analyses. For the case of WH — > ivbb, an isolated lepton 
i (electron or muon) and at least two jets are required, with one or more jets tagged as originating from 6-quarks. Two 
orthogonal tagging criteria are defined, one being an exclusive single-tag (ST) and the other a double-tag(DT) selection. 
These events must also display a significant imbalance of momentum in the plane transverse to the beam axis (referred 
to as missing energy or $ T ). Events with additional isolated leptons are vetoed. For the ZH — ► vvbb analyses, the 
selection is similar, except all events with isolated leptons are vetoed and stronger multijet background suppression 
techniques are applied. As there is a sizable amount of WH — > Ivbb signal that can mimic the ZH — > vvbb final state 
when the lepton is undetected, the D0 analyses include this as a separate search, referred to as WH — > ivbb. CDF 
includes this as part of the acceptance of the ZH — > vvbb search. In the WH — > Ivbb and ZH — > vvbb analyses, the 
final variable used for setting a cross section limit is the dijet invariant mass. The ZH — > £ + £~bb analyses require 
two isolated leptons and two jets, wherein the CDF analysis requires at least one of the jets to be ^-tagged while the 
D0 analysis uses only events with two 6-tags. For the D0 analysis the dijet invariant mass is used for setting limits 
while CDF uses the output of a 2-dimensional Neural Network to discriminate between signal and background. For 
the H — > W + W~ analyses, a large $ T and two opposite-signed, isolated leptons (electrons or muons) are selected, 
defining three final states (e + e~, e ± fi T , and The presence of neutrinos in the final state prevents the direct 

reconstruction of the Higgs mass, and the final variable used is the difference in (p between the two final-state leptons. 
The D0 experiment also contributes three WH — > WW + W~ analyses, where the associated W boson and the same- 
charged W boson from the Higgs decay semi-leptonically, thereby defining six final states containing all decays of the 
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* The Tevatron New-Phenomena and Higgs working group can be contacted at TEVNPHWG@fnal.gov. More information can be found 
at http://tevnphwg.fnal.gov/ 
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TABLE I: The luminosity, expected signal, expected background, and observed data for the CDF analyses. Also included is 
the Higgs mass for which each set of numbers are derived. For the ZH —* vvbb analysis the signal events are given separately 
for WH and ZH events. The numbers of expected events are determined for < rrijj <200 GeV/c 2 for H — > bb analyses, and 
< A(p(£ 1} £ 3 ) < 7T for H -> W+W~ analyses. 















WH -> £vbb 


ZH -> vvbb 


ZH ->£+rbb 


H -» W+W~ -> f ± v£^v 




DT(ST) 


DT(ST) 






Luminosity ( pb _1 ) 


1000 


1000 


1000 


360 


Expected Signal (evts) 


0.50 (1.47) 


0.25+0.22 (0.77+0.64) 


0.64 


0.84 


Expected Background (evts) 


40.2 (366.3) 


19.6 (309.9) 


103.7 


13.8 


Data (evts) 
m H (GeV/c 2 ) 


36 (390) 


24 (333) 


104 


18 


115 


115 


120 


160 


Reference 


[12] 


[13] 


[14] 


[15] 



TABLE II: The luminosity, expected signal, expected background, and observed data for the D0 H — > bb analyses. The number 
of expected signal, background, and data are determined for < rrijj <200 GeV/c 2 . 





WH -> evbb 
DT(ST) 


WH -> iMvbb 
DT(ST) 


WH -> Ivbb 
DT(ST) 


ZH -> vvbb 
DT(ST) 


ZH -» £+Tbb 


Luminosity ( pb _i ) 
Expected Signal (evt) 
Expected Background (evts) 
Data (evts) 
m H (GeV/c 2 ) 
Reference 


371 
0.19 (0.22) 
14.6 (58.3) 
14 (57) 

115 

M 


385 
0.11 (0.12) 
10.2 (46.8) 
8 (48) 

115 

[16] 


260 
0.18 (0.20) 
25.2 (87.3) 
23 (98) 

115 

[17] 


260 
0.24 (0.26) 
25.2 (87.3) 
23 (98) 

115 

U2 


320-389 
0.11 
11.0 
14 
115 
[18] 



third W boson (of opposite charge). In this case of this analysis, the final variable is a likelihood discriminant formed 
from several topological variables. 

All Higgs signals are simulated using PYTHIA v6.202[3(, using CTEQ5L[4] leading orderparton distribution func- 
tions. The signal cross sections are normalized to next-to-next-to-leading order calculations yjjllg], and branching ratios 
from HDECAYf?]. For both CDF and D0 , events from multijet (instrumental) backgrounds ("QCD production") 
are measured in data, but with different methods. For CDF, inherent backgrounds from other SM processes were 
generated using PYTHIA, ALPGEN]8], and HERWIG[9] programs. For D0 , inherent backgrounds were generated 
using PYTHIA, ALPGEN, and COMPHEP[l0|, with PYTHIA providing parton-showering and hadronization for 
all the generators. Background processes were normalized using either experimental data or next-to-leading order 
calculations from MCFMp]]. 

Values of integrated luminosity, and expected signal, expected background, and observed events are given in Ta- 
ble [T] for CDF analyses and Tables HIllIIII for D0 analyses. The numbers of events are integrated in the range 
< rrijj <200 GeV/c 2 for H — > bb analyses, < Acp(^i,£ 2 ) < tt for H — > W + W~ analyses, and < Discrim < 1 for 
the WH — > WW + W~ analyses. The tables also include the value of the Higgs mass for which each set of numbers is 
derived. 



TABLE III: The luminosity, expected signal, expected background, and observed data for the D0 WH — ► WW + W analyses. 
The number of expected signal, background, and data are determined for < Aip(£i,£2) < 7r. 



WW+W~ -> e ± ve ± v WW + W~ -> e ± v l i ± v WW + W~ -* ^v^i 

Luminosity ( pb" 1 ) 384 368 363 

Expected Signal (evts) 0.043 0.101 0.066 

Expected Background (evts) 15.4 7.0 12.5 

Data (evts) 15 7 12 

m H (GeV/c 2 ) 155 155 155 

Reference [19] [19] [19] 



3 



TABLE IV: The luminosity, expected signal, expected background, and observed data for the D0 H — > W + W analyses. The 
number of expected signal, background, and data are determined for < A(p(£i,£2) < 7T. 



H — > W + W~ — > e + ve~v H^W + W~^e ± v^u H — » W + W~ — > fi + ufi~v 

Luminosity ( pb" 1 ) 950 950 930 

Expected Signal (evts) 0.64 1.50 0.54 

Expected Background (evts) 11.4 28.1 10.5 

Data (evts) 11 18 10 

m H (GeV/c 2 ) 160 160 160 

Reference [20] [20] [21] 



III. COMBINATION PROCEDURES 



To gain confidence that the final result does not depend on the details of the statistical formulation, we combine 
all separate results using both a Bayesian and a Frequentist approach. In both methods, distributions in the final 
variables arc binned according to their experimental resolution, rather than as single integrated values. Systematic 
uncertainties enter as uncertainties on the expected number of signal and background events in each analysis. Both 
methods use likelihood calculations based upon Poisson probabilities. 



A. Frequentist Method 

The Frequentist technique relies on the CL S method, using a log-likelihood ratio (LLR) as test statisticflj], as given 
by: 



N 

LLR n = 2j2(s i -n l Log{l + s l /b l )) (1) 
»=i 

where n denotes the hypothesis being tested (e.g., background-only or observed data) and the sum runs over the 
number of bins (and/or analyses) being combined. The value of CL S is then defined as the normalization of the 
signal+background hypothesis (CL s+ t,) by the background-only hypothesis (CLj,). This construction reduces the 
ambiguity of "unphysical" results (e.g., negative cross section limits) and separates properties of the estimator from 
that of the quantity being probed. 



B. Bayesian Method 

Because there is no information on the production cross section for the Higgs, the Bayesian technique 0] assigns a 
flat prior for the total number of Higgs events. For a given Higgs mass, the combined likelihood is a product of the 
likelihoods in the individual channels, each of which is a product over histogram bins: 

Nc Nbins 

C(R, s, b\n) - J] J] (2) 

i=i j=i 

where the first product is over the number of channels (Nc), and the second product is over histogram bins containing 
riij observed events, either in dijet mass for WH and ZH, in Aip of two leptons for H — > W + W~, or in the 
likelihood discriminant for WH — > WW + W~ events. The parameters that contribute to the expected bin contents 
are [Aij = R x Sij + bij for the channel i and the histogram bin j. The posterior density function is then integrated 
over all parameters (including correlations) except for R, and a 95% credibility level upper limit on R is estimated by 
calculating the value of R that corresponds to 95% of the area of the resulting distribution. 
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TABLE V: The breakdown of systematic uncertainties for each individual CDF analysis. All positive-signed uncertainties 
within a group are considered 100% correlated across channels. Values with negative signs are considered uncorrelated. 



Source WH -> tvbb ST 


W# -> DT ZH 


-> z^&& ST 


ZH -> vvbb DT 


zh -> 


# -> iy + w- 


Luminosity (%) 


6.0 


6.0 


6.0 


6.0 


6.0 


6.0 


6-Tag Scale Factor (%) 


5.3 


16.0 


8.0 


16.0 


8.0 


n/a 


Lepton Identification (%) 


2.0 


2.0 


2.0 


2.0 


1.4 


3.0 


Jet Energy Scale (%) 


3.0 


3.0 


6.0 


(1.0-20.0) 


(1.6-20.0) 


1.0 


I(S)R+PDF (%) 


4.0 


10.0 


4.0 


5.0 


2.0 


5.0 


Trigger (%) 


0.0 


0.0 


3.0 


3.0 


0.0 


0.0 


Z + h.f. Shape (%) 


n/a 


n/a 


n/a 


n/a 


-20 


n/a 


Backgrounds 


W/Z+HF(I) (%) 


33.0 


34.0 


12.0 


12.0 


40.0 


n/a 


W+HF(II) (%) 








-10.0 


-42.0 





n/a 


Z+HF(II) (%) 








-6.0 


-19.0 





n/a 


Mistag (%) 


22.0 


15.0 


17.0 


17.0 


17.0 


n/a 


Top I (%) 


13.5 


20.0 


12.0 


12.0 


20.0 


n/a 


Top II (%) 


n/a 


n/a 


-2.0 


-3.0 


n/a 


n/a 


QCD (%) 


17.0 


20.0 


-10.0 


-44.0 


-50.0 


n/a 


Diboson I (%) 


16.0 


25.0 


12.0 


12.0 


20.0 


11.0 


Diboson I (%) 


n/a 


n/a 


-5.0 


-10.0 


n/a 


n/a 


Others (%) 


n/a 


n/a 


n/a 


n/a 


n/a 


-(12.0-18.0) 



C. Systematic Uncertainties 

Systematic uncertainties on background are generally several times larger than the expected signal, and are therefor 
important in the limit calculation. Each systematic uncertainty is folded into the signal and background expectations 
assuming Gaussian distributions. The Gaussian values are sampled once for each Poisson Monte Carlo (MC) trial 
(pseudo-experiment). Correlations among systematic sources are carried through in the calculation. Systematic 
uncertainties differ between experiments and analyses, both for signal and background sources. Detailed discussions 
of these issues can be found in the individual analysis notes [l|, H[ • Here we will consider only the largest contributions 
and correlations among and within the two experiments. 

1. Correlated Systematics 

The uncertainty on the measurement of the integrated luminosity is 6% (CDF) and 6.5% (D0 ). Of this value, 4% 
arises from the uncertainty on the inelastic pp scattering cross section, which is correlated between CDF and D0 . 
The uncertainty on the production rates for top-quark processes (ti and single-top) and electroweak processes (WW , 
WZ, and ZZ) are also taken as correlated between the two experiments. As the methods of measuring the multijet 
(QCD) backgrounds differ between CDF and D0 , there is no correlation assumed for this uncertainty. 

2. CDF Systematics 

The dominant systematic uncertainties for the CDF analyses are shown in Table [Vj For H — > bb, the largest 
uncertainties on signal arise from the 6-tagging scale factor (5.3-16%), jet energy scale (1-20%), and MC modeling 
(2-10%). For H — > W + W~ , the largest contributing uncertainty comes from MC modeling (5%). For inherent 
backgrounds, the uncertainties on the expected rates range from 11-40% (depending on background). Because the 
largest background contributions are measured using data, these uncertainties are treated as uncorrelated for the 
H — > bb channels. For the H — > W + W~ channel, the luminosity uncertainty is taken to be correlated between signal 
and background. The differences between treating the remaining uncertainties to be correlated or uncorrelated is less 
than 5%. 
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TABLE VI: List of leading correlated systematic uncertainties for the D0 analyses. The values for the systematic uncertainties 
are the same for the ZH — > vvbb and WH — > f.vbb channels. All uncertainties with a common origin are considered 100% 
correlated across channels. The correlated systematic uncertainty on the background cross section (a) is itself subdivided 
according to the different background processes in each analysis. 



Source 


WH -> evbb DT(ST) 


WH -> iivbb DT(ST) 


H -> W+W~, WH -> WW+W~ 


Luminosity (%) 


6.5 


6.5 


6.5 


Jet Energy Scale (%) 


4.0 


5.0 


3.0 


Jet ID (%) 


6.8 


6.8 





Electron ID (%) 


6.6 





2.3 


Muon ID (%) 





4.9 


7.7 


6-Jet Tagging (%) 


8.5(5.0) 


8.5(5.0) 





Background a (%) 


6.0-19.0 


6.0-19.0 


6.0-19.0 


Source 


ZH -> vvbb DT(ST) 


ZH -> e+e-bb 


ZH -» n + n~bb 


Luminosity (%) 


6.5 


6.5 


6.5 


Jet Energy Scale (%) 


6.0 


7.0 


2.0 


Jet ID (%) 


7.1 


7.0 


5.0 


Electron ID (%) 





8.0 





Muon ID (%) 








12.0 


&-Jet Tagging (%) 


9.6(6.7) 


12.0 


22.0 


Background a (%) 


6.0-19.0 


6.0-19.0 


6.0-19.0 



3. D0 Systematics 

The dominant systematic uncertainties for D0 analyses are shown in Table I VII The H — > bb analyses have 
an uncertainty on the ^-tagging rate of 5-7% per tagged jet. These analyses also have an uncertainty on the jet 
measurement and acceptance of 6-9% (jet identification or jet ID, energy scale, and jet smearing). For the H — > 
W + W~ and WH — ► WW + W~ , the largest uncertainties are associated with lepton measurement and acceptance. 
These values range from 3-6% depending on the final state. The largest contributing factor to all analyses is the 
uncertainty on the cross sections for inherent background at 6-15%. These systematics are assumed to apply to both 
signal and background processes. All systematic uncertainties arising from the same source are taken to be correlated 
between signal and background, as detailed in Table fVTl 

IV. COMBINED RESULTS 

Using the combination procedures outlined above, we extract limits on SM Higgs boson production a x BR* (H — > X) 
in pp collisions at ^/s = 1.96 TeV. Figure Q] displays the log- likelihood ratio distributions for the combined analyses 
as a function of ma. Included are the results for the background-only hypothesis (LLRf,), the signal -(-background 
hypothesis (LLR s+b ), and the observed data (LLR ofcs ). The shaded bands represent the 1 and 2 standard deviation 
(a) departures for LLR&. These distributions can be interpreted as follows: 

• The separation between LLRb and LLR s +b provides a measure of the overall power of the search. This is the 
ability of the analysis to discriminate between the s + b and b— only hypotheses. 

• The width of the LLR/, distribution (shown here as 1 and 2 standard deviation bands) provides an estimate of 
how sensitive the analysis is to a signal-like fluctuation in data, taking account of the presence of systematic 
uncertainties. For example, when a 1-er background fluctuation is large compared to the signal expectation, the 
analysis sensitivity is thereby limited. 

• The value of LLR D h s relative to LLR s+ f, and LLR^ indicates whether the data distribution appears to be more 
signal-like or background-like. As noted above, the significance of any departures of LLR Q fc s from LLR;, can be 
evaluated by the width of the LLRt, distribution. 

To facilitate model transparency and to accommodate analyses with different degrees of sensitivity, we present our 
results in terms of the ratio of limits set to the SM cross sections as a function of Higgs mass. A value of < 1 would 
indicate a Higgs mass excluded at 95% CL. The expected and observed 95% upper limit ratios to the SM cross section 



6 



for the combined CDF and D0 analyses are shown in Figure [2j The observed and expected limit ratios are listed for 
selected Higgs masses in Table IViTI with observed(expected) values of 10.4(7.6) at m# = 115 GeV/c 2 and 3.8(5.0) at 
m H = 160 GeV/c 2 . 

These results represent an improvement in search sensitivity over those obtained for individual experiments, which 
have set observed(expected) limit ratios of 16.3(16.7) for D0 and 12.8(9.1) for CDF at m# = 115 GeV/c 2 and of 
4.3(5.9) for D0 and 11.4(10.2) for CDF at m H = 160 GeV/c 2 . 
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TABLE VII: Expected and observed 95% CL cross section ratios for the combined CDF and D0 analyses. 



100 GeV/c 2 


115 GeV/c 2 


120 GeV/c 2 


140 GeV/c 2 


160 GeV/c 2 


180 GeV/c 2 


200 GeV/c 2 


Expected 6.2 


7.6 


8.7 


9.3 


5.0 


7.5 


15.5 


Observed 8.5 


10.4 


11.1 


8.8 


3.8 


6.1 


12.3 




FIG. 1: Log-likelihood ratio distribution for the combined CDF and D0 analyses. Shown in the plot are the LLRt (background- 
only hypothesis), LLR S+ ;, (signal+background hypothesis), LLR (, S (observed LLR value), and the 1-a and 2-a bands for the 
LLRj, distribution. 
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FIG. 2: Expected and observed 95% CL cross section ratios for the combined CDF and D0 analyses. Also shown are the 
expected 95% CL ratios for the CDF and D0 experiments alone. 



